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abstract 



Large-A^ QCD with heavy adjoint fermions emulates pure Yang-Mills theory at long distances. 
We study this theory on a four- and three-torus, and analytically argue the existence of a large- 
small volume equivalence. For any finite mass, the center-symmetry unbroken phase exists at 
sufficiently small volume and this phase can be used to study the large volume limit through 
the Eguchi-Kawai equivalence. A finite-temperature version of volume independence implies that 
thermodynamics on x can be studied via a unitary matrix quantum mechanics on S^, by 
varying the temperature. To confirm this nonperturbatively, we numerically study both zero- and 
one-dimensional theories by using Monte-Carlo simulation. The order of finite- corrections turns 
out to be 1/A^. We introduce various twisted versions of the reduced QCD which systematically 
suppress finite- A corrections. Using a twisted model, we observe the confinement/deconfinement 
transition on a 1'^ x 2 lattice. The result agrees with large volume simulations of Yang-Mills 
theory. We also comment that the twisted model can serve as a nonperturbative formulation of 
the noncommutative Yang-Mills theory. 



E-mail address : aze@gauge.scphys.kyoto-u.ac.jp 
^ E-mail address : masanori.hanada@weizmann.ac.il 
E-mail address : unsal@slac.stanford.edu 
E-mail address : ran.yacoby@weizmann.ac.il 



Contents 

1 Introduction 

1.1 Results 

2 Review of original Eguchi-Kawai proposal 

2.1 Perturbative fluctuations and one- loop potential 

2.2 Nonperturbative fluctuations and size of eigenvalue bunch . . 

3 Eguchi-Kawai equivalence for massive QCD(Adj) 

3.1 QCD(Adj) at zero temperature on at small L 

3.1.1 iVf = (bosonic) 

3.1.2 iV^ = 1/2 (single Maj orana ) 

3.1.3 Nj^ > 1: Uniformization of eigenvalue distribution . . 

3.2 QCD(Adj) at finite temperature on asymmetric T'^ x 5^ . . . 

3.2.1 Nj? = (bosonic) 

3.2.2 Nf = 1/2 (single Majorana) 

3.2.3 Nf >1 

4 Lattice model and Monte-Carlo simulation 

4.1 Single-site theories 

4.1.1 QCD(Adj) on 1^ lattice 

4.1.2 TQCD(Adj) at zero temperature 

4.1.3 Numerical results for Nj' = 1 

4.2 One-dimensional lattice: Reduced model at finite temperature 

4.2.1 (T)QCD(Adj) at finite temperature 

4.2.2 Monte Carlo simulation for NjP = 1 

5 Stabilizing noncommutative Yang-Mills theory 

6 Conclusions and Discussions 



A QCD(Adj) on continuum at finite-L and one-loop analysis 

A.l IR singularities, conditional convergence, and noncommutative saddles 



1 Introduction 



Recently Yang-Mills (YM) theory with adjoint fermions, QCD(Adj), has attracted much in- 
terest. The main impetus behind this is a network of exact large-A^ equivalences. The starting 
point is the large- orientifold equivalence, which states that the bosonic subsector of this theory 
is equivalent to the charge-conjugation even subsector of QCD with fermions in antisymmetric 
representation [QCD(AS)] [1], provided symmetries defining the neutral subsectors are not spon- 
taneously broken [2]. QCD(AS) reduces to the ordinary QCD with fundamental quarks when 
N = 3, and is a natural large- A^ generalization thereof |3]0 The second important link is an orb- 
ifold equivalence: when quarks are massless or light with respect to strong scale Aqcd, QCD(Adj) 
compactified on a Euclidean four-torus exhibits volume independence, thanks to its unbroken 
prerequisite (center and translational) symmetries at any radius [7]. Thus, through the Eguchi- 
Kawai (EK) equivalence [8], one can study large-A^ QCD on R"^ by using a unitary matrix model 
on a single-site lattice. QCD(Adj) also provides new insights into gauge dynamics, especially on 
small X . This theory exhibits new nonperturbative phenomena, most strikingly confinement 
due to magnetic bions, a new class of non-self-dual topological excitations [Hi [10], distinct from 
monopoles and instantons. 

The statement of the volume independence [TJ [SJ [TU [T^j is as follows. Consider SU(N) gauge 
theories on M^, with or without fermions, toroidally compactified on a four manifold M.^~'^ x T'^. 
For simplicity, we assume the matter fields are in adjoint representation and, hence, the theory 
has a global {'^nY center symmetry, described most easily as gauge rotations aperiodic up to an 
element of the center group. The order parameters of this symmetry are Wilson lines wrapping 
distinct toroidal cycles. The observables singlet under the center transformation constitute the 
neutral sector. The volume independence states that dynamics of the neutral sector observables 
is independent of the size of the torus provided the center symmetry and translational invariance 
are not spontaneously broken. Among such observables are nonperturbative mass spectrum, free 
energy densities, and deconfinement transition temperature, just to count a few. This is clearly an 
extraordinarily well justified reason to study aspects of the small volume, large- A^ gauge theories. 

In fact, center symmetry is spontaneously broken in most examples. In the original Eguchi- 
Kawai model [8], which is a one-point reduction of Wilson's bosonic lattice gauge theory, the break- 
down can be shown by one-loop calculation around a diagonal background [IB]. The quenched 
[131 and twisted [15] modifications of the Eguchi-Kawai model were proposed to preserve the 
symmetry, but after two decades, it is now understood that both modifications fail nontrivially 
due to nonperturbative effects [E], [13 EHl 19] | 1 

Similarly, in any gauge theory compactified 
thermally on x S"^, or on a torus with at least one thermal boundary condition, center symme- 
try breaks spontaneously in the high temperature deconfined phase, and volume independence is 
only valid in the low temperature confined phase, above a critical volume jl2j . 

Kovtun, one of us (M.U.), and Yaffe, motivated by the quantitative differences between ther- 

^The phenomenology of QCD(AS) is examined in [4l[5l[6]. 

^The failure of these modifications can be cured by introducing supersymmetry [191120] , Reference [21] proposed 
a concrete way to construct 4d A/" = 4 super Yang-Mills theory by using the Eguchi-Kawai equivalence, which 
preserves 16 supersymmetries. Reference [22] studied the Eguchi-Kawai reduction in the strong coupling domain of 
A/" = 4 super Yang-Mils by using AdS/CFT and D-branes. 

''Recently, a new limiting procedure for the twisted Eguchi-Kawai model, which aims to prevent center breaking, 
has been proposed |23j . 
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mal and circle (nonthermal) compactifications, showed that if light or massless adjoint fermions 
endowed with periodic boundary conditions are added to Yang-Mills theory, then the center sym- 
metry is stabilized at small volume dynamically [7]. (Also see the discussion in Refs.|24l I25[ I26j.) 

For heavy fermions, the infrared physics of QCD(Adj) on emulates the bosonic Yang-Mills 
theory. Since heavy fermions are also capable of restoring center symmetry at sufficiently small 
volume, this may provide an opportunity for a working Eguchi-Kawai reduction for an "almost" 
Yang-Mills theory. However, this is not straightforward. When one dimension is compactified, 
QCD(Adj) with massive fermions on small x M'^ exhibits an intricate phase structure. This 
can be deduced from a one- loop effective action of the Wilson line \25\ I27j . simulations on an 
asymmetric torus mimicking x [25] , and studies on x which also mimics x 

due to topological reasons, as explained in [30]. In all these cases, the Z^r symmetry along is 
intact at large radius, and as one decreases the radius, it breaks down completely at some critical 
point and then gradually restores to various subgroups of the center symmetry. Zjy symmetry is 
restored fully only at mLN ~ few, where L is the compactification radius and m is the fermion 
mass in continuum. This is troubling because the volume independence (strong coupling, non- 
Abelian confinement) domain is LA^Aym ^ 1 [21]) whereas, for heavy fermions, m > Aymj 
the first condition implies LA^Aym ^ few which is a volume dependent, weak-coupling Abelian 
confinement domain [31j . 

One might expect a similar pattern for a single-site lattice model based on this weak-coupling 
intuition. However, recent important work of Bringoltz and Sharpe shows that (Z^r)^ remains 
intact in a rather generous "funnel," in the fermion mass, lattice coupling plane, covering the 
continuum limit of Yang- Mills theory [32] (also, see [33]), corresponding to the limit where the 
bare fermion mass is larger than cutoff scale. In particular, they observe a (Z^r)^ restoration 
at ma ~ 0{N^), where m is bare lattice mass and a is lattice spacing. Why and how the full 
center-symmetry restoration takes place at ma ~ 0{N^) is the main theoretical problem that we 
wish to address in this work. 



1.1 Results 

To set the notation, we first express the action of continuum QCD(Adj) on a four manifold: 



j^t I d^xTr[^F'^, + J2i;f{p + m)i;f], (1) 



where tpf are Dirac fermions with mass m. (Generalization to different values of masses is straight- 
forward.) In general, we will consider this theory on four-torus x 5^, with sizes L and /3, 
respectively. If we impose the antiperiodic boundary condition on fermions along the temporal 
direction, this action describes the finite-temperature system and /? corresponds to the inverse 
temperature. For periodic boundary conditions on fermions in all directions, we consider (3 = L, 
a symmetric four-torus. 

We show that, with periodic boundary conditions, {^n)'^ symmetry is not broken at sufficiently 
small-L, although it can be broken at some intermediate volume. (More precisely, (Z^r)^ symmetric 
and broken phases can coexist, while quantum tunneling between them is suppressed in the large- 
N limit.) The argument, which will be quantified in Sec l3.1l is very simple - although the one-loop 
effective action suggests the existence of the attraction between Wilson line eigenvalues at small 
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separation, the eigenvalues spread due to nonperturbative quantum fluctuations [l2], and the one- 
loop calculation can be trusted only at large separations where it leads to repulsion. The estimates 
of nonperturbative fluctuations are outside the reach of one-loop perturbation theory and often 
overwhelm the implications of one-loop analysis. Because of nonperturbative effects, we find that 
the full center restoration takes place at mL ~ 0{N^), which is compatible with LA^Aym ^ 1 for 
m > Aym- 

Our main results are 

1. Small-large volume equivalence: The QCD(Adj) with heavy fermions of mass m on M^, 
or with radii L > A~^, is equivalent to the theory on a small T^, with radii mL < 0{N^), 
due to unbroken (Z^r)^ center symmetry associated with the cycles on T^. This is a small- 
large volume equivalence with an intermediate center-symmetry broken phase, where volume 
independence is not valid. This is depicted in the following figure. 

■ > • 

mL ~ 1 AL ~ 1 00 



2. Finite-temperature equivalence: The theory on x at finite temperature is equiva- 
lent to the one on x provided the {^nY center symmetry associated with the cycles on 

is not spontaneously broken. The phase transition in the thermodynamic x limit 
can be studied by using a large- reduced model on small x by dialing /3. This form of 
equivalence is also useful for Hamiltonian formulation and extraction of the nonperturbative 
spectrum of the theory. 

3. Twisted QCD: In compactified QCD(Adj), finite-A^ corrections turn out to be order 1/A^, 
as opposed to perturbative expectation on [3l], which is order l/N"^. There are two 
plausibly related explanations for this behavior, whose footprints can be seen in finite- volume 
perturbation theory. In a weak-coupling center-symmetric background, the volume is only 
enhanced by a factor of A^, and the effective volume is V^g ~ NV . Finite-A^ corrections 
should scale as finite-volume corrections. The other is, in compact space, one cannot gauge 
away zero-momentum modes. Typically, there are order N bosonic and fermionic zero modes, 
which may generate nonperturbative 1/A^ effects j35| . Both problems can simultaneously 
be cured and 1/N corrections can systematically be improved by using twisted boundary 
conditionsEI We refer the latter as TQCD(Adj) as per [15]. 

Our results have interesting spin-offs for noncommutative theories, phase transition in pure 
Yang-Mills theory and orientifold equivalence. Adding massive or massless adjoint fermions to the 
twisted Eguchi-Kawai model (TEK) model cures the global instability of the model \n \ [T9 t [TH]. 
Therefore, our formulation can be used to provide a nonperturbative definition of noncommuta- 
tive bosonic Yang-Mills theory and noncommutative QCD(Adj). By using the reduced matrix 

*The relation between "effective" volume and is discussed in many places, for a review, see [3^. The obser- 
vations that twisted boundary conditions can be used to systematically reduce finite- volume corrections is given 
in Ref.[37] and ii) to lift bosonic and fermionic zero modes is given in Ref . [38] . We give analytic and numerical 
evidence suggesting that the two are indeed related. 
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model for TQCD(Adj) with very heavy adjoint fermions on 1'^ x 2 lattices, we observe the con- 
finement /deconfinement transition at bare coupling be = 0.32 — 0.33. Large volume simulations 
for pure Yang-Mills theory give similar results, in accordance with the finite-temperature version 
of equivalence (For example, SU (8) YM theory simulated on 10^ x 5 lattices gives, in the extrapo- 
lated infinite volume limit, be ~ 0.34 [39] (see also [10]). The small difference is due to 1/A'' effects, 
the existence of heavy fermions in our reduced model and the difference of numbers of sites along 
the temporal direction). We hope that the unitary matrix model can be used to gain insight into 
the nature of the deconfinement transition of the infinite volume theory. 

Finally, combined with the large- orientifold equivalence between (T)QCD(Adj) and QCD(AS), 
thermal properties of the QCD(AS) can also be studied by using the large- A^ reduced model for 
TQCD(Adj) at finite temperature P [3 [30] Jl Our results for massive QCD(Adj) can be regarded 
as a first step towards this direction of research. 

2 Review of original Eguchi-Kawai proposal 

In this section, we give a brief review of the Eguchi-Kawai equivalence. Our aim is not to repeat 
the original works on the subject; rather we wish to point out two important ingredients which 
will be useful. These are perturbative and nonperturbative quantum fluctuations in the matrix 
model. 

The EK model is the dimensional reduction of Wilson's lattice gauge theory action down to a 
1^ lattice, i.e., a one-site matrix model. The action of the reduced model is 

Sod = -2bNRe Tr ( V^V.V^Vj ] , (2) 

where {fi = 1, 2, 3, 4) are unitary matrices, b is the bare inverse 't Hooft coupling constant. To 
probe continuum physics, b should be chosen appropriately depending on the lattice spacing a. 
This action has a (Z^r)^ global center symmetry 

y^^^2.^n,/Ny^ (^^ = Q, 1, • • • , A - 1). (3) 

The reduced model ([2]) is equivalent to the translationally invariant subsector of the lattice theory 
with an arbitrary number of sites, provided center and translational symmetries are not broken. 
These are the necessary and sufficient symmetry realization conditions for the large-A volume 
independence. 

2.1 Perturbative fluctuations and one-loop potential 

The realization of center symmetry can be determined by integrating out weakly coupled pertur- 
bative modes in the background of diagonal, commuting Wilson lines: 

y^ = diag(e<,--- ,e<), [V^,V,]=0. (4) 

'''Earlier work on deconfinement transition in large-A'^ QCD(Adj) and QCD(AS) in the weak-coupling limit of 
small S*'^ X 5*^ showed that the critical temperatures agree [30j . up to corrections. 
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The resulting one-loop action reads (we parallel the discussion in the one-site matrix model and 
continuum theory on small T^) 



S 



1— loop 



sm 



^ 2^a<b l^kez* iog 12^^=1 7^ 



on 1 — lattice 



on — continuum 



(5) 



fi Note that one-site_ matrix model and have the same symmetry properties, 



where Of!' 

both invariant under 0^ — ?• 0^ -|- 27r|f| It is also convenient to rewrite ([5]) by using the Wilson 
lines. This can be done by a Fourier transformation for the former and by a Poisson resummation 
for the latter: 



'S'l-loop[^l, ■ ■ ■ , V4] — < 



*\{0} 



n4\|2 



N 



on 1^ — lattice 



on T'* — continuum 



(6) 



where 



P{n) 



I 



da 



a 



e-«"I„,(2a)...I„,(2«) 



(7) 



and In{2a) is the modified Bessel function of the second kind; InC^a) = (27r)~^ 
In the one-site model, Fourier coefficients are slightly more complicated, but the main result is 
the same as in continuum T^. The basic point is the negative definiteness of the integral in ((71). 
Moreover, at a large-winding number, the integral can be evaluated analytically by localization, 
and converges to the continuum result quickly: 



P{n) < 0, Vn G \ {0}; 



P{n) 



1 1 



for Inl ^ 1 , 



(8) 



The one-loop action ([5|) has IR singularities whenever two (or more) eigenvalues coincide. This 
is also manifest in ([6]); both series are conditionally convergent at a large- winding number, and 
exhibit logarithmic divergence for coinciding eigenvalues. These issues are discussed throughly in 
the Appendix |X1 Physically, at these points, there are extra massless degrees of freedom which 
should not have been integrated out. In other words, the zeroth order assumption that one can 
expand the fluctuations around commuting saddles (Uj) is not always correct. As discussed in the 
Appendix |Al the IR divergence can be regularized and meaningful results can be extracted from 
([5]). The result is, of course, well-known; ([5]) generates an eigenvalue attraction, or in Q the 
"masses" for Wilson lines are all negative and, consequently, eigenvalues clump. However, neither 
implies that all the eigenvalues are coincident. Because of nonperturbative effects, the eigenvalues 
spread. In order to estimate the size of the eigenvalue bunch, we may study the theory around 
one of its global minima, ~ 1. 



The relation between the continuum one-loop effective action on T" and the one-loop effective action for one-site 
theory given in Eq. ^ is analogous of the one between the action of the XY model and its Villain form. 
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2.2 Nonperturbative fluctuations and size of eigenvalue bunch 

Parametrizing w g«a^M ^ observe that around the minimum of the one- loop potential, all the 
fluctuations are quartic (as opposed to being quadratic). The action expanded around one of the 
minima, the ~ 1 configuration gives 

Sod = ^rrf-i[X^Xn'Y (9) 




4^ 



where 



1 



on 1^ — lattice 



^od=<| (10) 
'^r4^ ~ — , \ ou — continuum 

is the zero-dimensional 't Hooft coupling. In the lattice, b is, of course, a parameter that one 
can choose at will. However, to probe continuum physics, it needs to scale as 6 ~ log(l/aA) by 
asymptotic freedom. 

Since the Hermitian matrix model is given in terms of noncompact matrices, it is not a priori 
guaranteed that the theory actually exists quantum mechanically. Reference [H] shows that the 
partition function of the theory defined through ([9]) does not exist for SU(2), but exists for SU(N) 
for > 3. In the next section, we will indeed see closely related Hermitian matrix models 
which do not exist for any N. Despite this, the Hermitian matrix model defined through ([9]) is 
useful. It can be employed to determine the interaction between eigenvalues and to estimate the 
root-mean-square fluctuations of X^^ matrices. (We also refer to this as the size of the eigenvalue 
bunch.) 

Scalar eigenvalues in Od theory are related to the phases 9 of the Wilson lines winding on 
temporal and spatial directions by 

1x1^ = 9^ (a = l,--- ,iV). (11) 

The 't Hooft coupling Aod has the dimension of (mass)^ and its value sets the typical mass scale 
of the Od theory In particular, typical fluctuation of the eigenvalues of the dynamical fields is 
set by this scale. Because of the generic noncommutativity of matrices , the relative positions 
of the eigenvalues make sense only when their separation is of order or larger than Aq^^ [12] . 

The one-loop effective action of the matrix model around the diagonal background can be 
calculated if all eigenvalue pairs are well-separated \xa — xi,\ ^ ^^od^^ corresponding to the case 
where all eigenvalues are weakly coupled. Integrating out massive 'W bosons" (off-diagonal 
elements) by using the background field gauge yields 

Si-loop[x'^b] = 2^1og|Xa - Xfop. (12) 

a<b 

Clearly, this is nothing but Eq.Q restricted to its Kaluza-Klein (KK) zero mode. Thus, as in the 
one-loop potential ([5]) of the full theory, the Hermitian model as well predicts eigenvalue attraction 
and also exhibits the same IR singularity whenever two eigenvalues coincide. 
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More importantly, the Hermitian matrix model ([9]) allows the determination of the typical size 
of eigenvalue fluctuations. In 't Hooft's large-A^ limit, the root-mean-square fluctuations AX are 



This result has various interesting implications. 

The analysis of continuum and weak-coupling (6—7-00) domain of lattice gauge theory 
exhibits similar behavior. Let L = a. Then, AX ~ 1/ [log(l/LA)]^^^ L < 1/L. At asymptotically 
small L, the size of the eigenvalue bunch is much smaller than the size 2tt/L of the dual torus where 
eigenvalues are residing. In this domain, clearly, nonperturbative fluctuations cannot overwhelm 
broken center symmetry. 

It is well-known that Eguchi-Kawai reduction holds in the strong coupling domain of lattice 
gauge theory, for < 6 < 6c = 0.19. This is a lattice domain unrelated to continuum physics, 
whereas the above analysis is valid for sufficiently large 6. What happens to AX if we decrease 
b to being of order few? Of course, this is unjustified, but demonstrates the trend of the root- 
mean-square fluctuations. Naive use of the AX ~ 1/6^/^0 shows that when 6^/^ is order one, the 
quantum fluctuations may be as large as the size of the dual torus. In this domain, eigenvalues 
feel the size of the compact space they live in and we expect the strong fluctuations to lead to 
center restoration. This is indeed the case. 

In this paper, we will study Yang-Mills theory with massive adjoint fermions with mass m. 
When m = 00, the discussion is the same as the original Eguchi-Kawai reduction. For ma ~ few, 
the perturbative-loop analysis indicates that only a finite subgroup of center symmetry would 
restore; however, such analysis does not take into account nonperturbative fluctuations. We 
propose that these fluctuations tend to uniformize the eigenvalue distribution rather quickly, even 
when ma ~ few, providing the resolution of the problem stated in the introduction. 

3 Eguchi-Kawai equivalence for massive QCD(Adj) 

In this section we explain why it is natural to expect that the Eguchi-Kawai reduction holds for 
the QCD(Adj) with massive adjoint fermions at sufficiently small volume. Our analysis is based 
on the one-loop potential for diagonal components of fields and estimates of nonperturbative 
quantum fluctuations. We start with the zero temperature case and then generalize to the finite 
temperature. 

3.1 QCD(Adj) at zero temperature on at small L 

Let us consider the QCD(Adj) at zero temperature on a symmetric four-torus, with size L. 
At sufficiently small L, gauge coupling at the scale of the compactification is small and we may 
analytically compute the one- loop effective action on T^. The reason for doing this computation, 
apart from trying to determine the center-symmetry realization, is twofold. One is we would like 
to compare with the one-site theory, given in ()33p . More importantly, we give evidence that some 
(not all) implications of one-loop action are, in full theory, overwhelmed by large nonperturbative 
quantum fluctuations, and therefore, incorrect. 




(13) 



9 



The one-loop effective action, in the Wilson line background induced by gauge and 
fermionic fluctuations with periodic boundary conditions can be written as 




5i-ioop[^,J = 2^ 2^ (2 log j-mj log I 2^ + m 

where Of = 91^-6^^, and m is the fermion mass. Similar to original EK model ([5]), this expression 
has IR singularities whenever two (or more) eigenvalues coincide, which we discuss throughly 
in the Appendix |XJ The theory may have different saddles which are expressed in terms of 
noncommuting matrices. The classiflcation of the saddles, using the techniques of Ref.[SB], of 
QCD(Adj) as a function of mass of the fermion is given in the Appendix [XI 

When we consider a regime where \e'^^\ <^ 2tt, we can split ()14p into the zero and non-zero- 
momentum contribution. The eigenvalue dynamics is dominated by the interactions between 
nearby eigenvalues and the effect of high KK modes is negligible. Therefore, we will use the trun- 
cated Hermitian matrix model (jlSp to gain an understanding of the typical eigenvalue fluctuations. 
Strictly speaking, the truncation can be justifled only when the (Zat)^ center symmetry is com- 
pletely broken and consequently there exists a clear separation of scales between the KK modes 
and zero modes. In a center-symmetric background, this is not the case. However, at large- A^, the 
states obtained by quantizing the theory on a center-symmetric vacuum fill the [0,2tt/L] energy 
range uniformly. If we consider a finite but small range \6"'^\ <C 27r, there are stifl 0{N^) states 
(in perturbation theory) in this interval. We may therefore use the Hermitian model to probe 
the interaction of nearby eigenvalues, and their fluctuations. The eigenvalue dynamics of the full 
theory is mimicked rather accurately by the truncated Hermitian matrix model. 

The truncation of the KK modes in ([T]) yields the zero-dimensional Hermitian matrix model 

TV" 

A^ 



Sod — T 



Tr ^[A^ + ^5/ (7M[A^ V/] + rm^f) j , (15) 



where Xod = \Ad{l/L)/L^ ~ 1/ log(l/LA)L^ is the zero-dimensional 't Hooft coupling, and A is 
the strong scale of QCD(Adj). 

The and m are the two scales in the Hermitian matrix model ()15p . In particular, for the 
m = oo theory, which corresponds to the original EK model discussed in SeciH the root-mean- 
square fluctuations of matrices are set by this scale AA(m = oo) ~ '^Jd ' With standard 't Hooft 
scaling, this is 0{N^) in the large-A limit. This is a crucial observation that will be important 
below. Because of the generic noncommutativity of A^j matrices , the relative positions of the 
eigenvalues make sense only when their separation is of order or larger than Aq^^ [42j . 

The one-loop effective action of the matrix model around the diagonal background can be 
calculated if all eigenvalue pairs are well-separated \xa — xi,\ ^ -^od^' corresponding to the case 
where all eigenvalues are weakly coupled. Integrating out massive "W bosons" (off-diagonal 



^It should be understood that the Hermitian model p5|) is used as auxiliary. In the strict L — Q limit, and for 
massless fermions, the partition function of the Hermitian model diverges. One-loop potential is unbounded from 
below as eigenvalue separations tend to infinity. At finite-L, the target space of eigenvalues is bounded, cannot 
run-off to infinity and the partition function of the finite-i theory exists. 
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elements) by using the background field gauge yields 

5i-ioop[<J = 2Y,iog\xa-Xb\^ -iNfY,log{\xa-Sb\^ + m^) . (16) 

a<b a<b 

Note that this is nothing but Eq. (|14p restricted to its KK zero mode, as expected. 

The effect of adjoint fermions, which is most manifest in (jl6p . is to generate repulsion among 
eigenvalues. At m = oo, since the fermions do not contribute, they have no impact on AX. 
For finite fermions' mass, the effect of fermions makes AX larger than that of the m = oo case. 
Although we will not prove this statement, it is easy to understand it on physical grounds. When 
fermion mass is zero, then the one-loop potential is unbounded from below only at large eigenvalue 
separation; hence AX{m = 0) = 00. From now on, we assume AX{m) > AX{m = 00) = X^^f. 
Somewhat conservatively, we use the smaller of the fluctuations in what follows. 

Let us now discuss the realization of (Z^v)^ symmetry based on the effective action (jl6p for 
the Hermitian matrix model. 

3.1.1 Nj? = (bosonic) 

The Nf = or m = 00 limits are same and are discussed in Secl2j We repeat the main 
result for convenience. The one-loop action leads to the mutual attraction of eigenvalues at large 
eigenvalue separation and hence, eigenvalues must clump. This implies broken center symmetry. 
However at small eigenvalue separation, one-loop approximation is not valid; the background of the 
commuting Wilson line saddles breaks down. Because of nonperturbative quantum fluctuations, 
the eigenvalues do not collapse to a point; rather the eigenvalue clump has a finite extent of the 
order AX ~ aJ^^ = [A4d(l/L)]i/^/L. Note that, although the size of the eigenvalue clump is 
suppressed with respect to 1/L at small A4rf(l/L) and hence center symmetry is broken, it scales 
as 0{N^) in the large- limit. This will be crucial later. 

3.1.2 Nj? = 1/2 (single Majorana) 

When Nj^ = 1/2 and m = 0, the theory is 4d A/" = 1 pure super Yang-Mills theory. In this 
case, the one- loop effective action (|16p vanishes and in fact, this is true to all loop orders due 
to super symmetry. Taking nonperturbative fractional instanton effects into account, the center 
is unbroken on x 5^ as discussed in [7]. In supersymmetric theories with supersymmetry 
preserving boundary conditions, it is expected that there are no phase transitions as the volume is 
varied [38]. At large- A^, the absence of phase transitions transmutes to exact volume independence 
[7] and unbroken (Ztv)^ center symmetry. However, it is also possible to construct a metastable 
center-broken sector [33], which becomes stable at large- A^. 

When m is fixed and nonzero, by taking L — )• 0, p6p is positive, and therefore, (^n)^ symmetry 
breaks down. If L is fixed and m — ?• 0, a center-symmetry preserving background exists for a finite 
range of m. This aspect is discussed throughly in Ref . |44j . This noncommutativity of limits 
requires care in drawing conclusions about this case. 

3.1.3 A^j* > 1: Uniformization of eigenvalue distribution 

The effective action (fT6|) predicts attraction at short distance |Ax| < m, and repulsion at 
long distance |Ax| > m. Then one may naively conclude that all eigenvalues clump and (Z^v)^ 
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symmetry is broken. However, one should notice that this effective action is vahd only at |A2;| > 
^Jd ■ When Aq^"^ < m attractive force emerges at Aq^'^ < |Ax| < m and hence (Zjv)^ is broken 
(left of FigdJ. When Aq^^ > m, however, it predicts only repulsion; eigenvalue fluctuation is of 

order Aq^^ and hence they cannot clump to small region where (jl6p predicts attraction (right of 
FiglU). 

At sufficiently small, but 0{N^) compactification radii L, we can always guarantee that non- 
perturbative quantum fluctuations overwhelm fermion mass, i.e., Aq^"^ = [A4d(l/L)]^/^/L > m. 
In this case, fermion mass is negligible and for the purpose of center-symmetry realization, the 
theory cannot be distinguished from the massless theory, for which center symmetry is unbroken. 
At such values of L, since the target space of eigenvalues is compact four-torus with size j^, 
the repulsion implies that the eigenvalues will uniformly distribute over T^. 

This is the sense in which lower dimensional nonperturbative quantum fluctuations help 
restoration of center symmetry at niL ~ few, as opposed to quantum field theory on x 
where full center restoration requires inLN ~ few. 




Figure 1: The scales in the problem. Left panel: m > Aq^ and mass is important. Right panel: 
m < Aq^^ and mass is negligible with respect to quantum fluctuations. We can always realize the 
latter case by taking the size of the four-torus sufficiently small, but still 0{N^). See the text for 
explanations. 

As we will see in the Appendix lAl the (Zat)^ broken phase which consists of k bunches of 
eigenvalues may exist when (Aod/^)^^^ f^i- However, the radius of each bunch is of order 
{Xod/kY^^ 1 which is smaller than the eigenvalue fluctuation in the (Z^v)^ symmetric phase (Aod)^^^, 
and hence even if the (Z^r)^ broken phase exists, tunneling to such a state is highly suppressed. 

The above argument explains the reason for the working of Eguchi-Kawai reduction in Refs.[32t 
[33] for heavy fermions. In [32], initial configuration for simulation corresponds to X^^ = (single 
bunch). It is unstable and collapses to (Z^r)^ symmetric phase. 

3.2 QCD(Adj) at finite temperature on asymmetric x 

We can generalize the argument of the previous subsection to the finite-temperature QCD(Adj) 
on an asymmetric four-torus x 5^. Circumferences of temporal and three spatial circles are 
taken to be (3 and L. On fermions, we impose antiperiodic and periodic boundary conditions 
along temporal and spatial circles, respectively. 
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Completely analogous to the discussion of Sec l3.1l at small L, the eigenvalue dynamics of the 
full theory is mimicked rather accurately by the truncated Hermitian matrix quantum mechanics 
with actionlf] 




N fP 1 1 \ 

Sid = x^J^ ^^'^n 2 (^*^')' - 4 [^'' ^']' + ^'^ (T^AV'a + ^i[X\iJa] + mV'a) j ,(17) 



where 

Ai. = ^. (18) 

The 't Hooft coupling Aid has the dimension of (mass)^ and its value sets the typical mass scale 
of the theory. In particular, typical fluctuation of the eigenvalues of the dynamical fields is given 
by this scale. The mapping between the scalar eigenvalues and the phases 9 of the Wilson loops 
is the same as in (jlip. but now only running over spatial directions. Since the S^ direction is 
compact, unlike its decompactification limit S^ — ?■ M, the Wilson line along that direction cannot 
be gauged away. We parametrize Wt = diag{e'^°'^^ , • • • , e*"^^), or At = log Wt- 

The one-loop effective action of the matrix quantum mechanics around the static diagonal 
background in the background field gauge can be computed in the weak-coupling domain, 



,1/3 

\d 



resulting in 



<C 1 or < 1 , (19) 



5'id,i-ioop(a,^) = 2^ ^ logii—- + {aa-ab)j + \xa - Xbl"^ \ 

a<b n=-oo \ V / / 

.ArDY- ^ 1 //27r(?i + l/2) ^ y _ ^,2 

a<bn=~oo \ V / 



(20) 



where and Xa represent diagonal components of the gauge field At and three scalars Xi, 
respectively. By subtracting a constant factor, (j20|) can be rewritten as [IS] (for reference, we 
quote both periodic and antiperiodic boundary conditions along /3 circle) 

'S'u.l-loop = 2^1og(cOsh(/3|fa - ffol) - COs(/3(aa - Ob))) 
a<b 

-ANf log (cosh [p^\Xa-Sb\'^+m?^ ^ COs(/3(Qa - Ob))) . (21) 

a<b 

From the expression of S^^ noop' oris might naively conclude that eigenvalues Xa and Ua coincide, 
because the one-loop action is negative infinity at that point. However, the analysis is valid only 
when condition ()19p is satisfied; thus drawing reliable conclusions requires care. 



®We repeat the same cautionary note as in Section [3TT] The classical theory written in Ea. (|17|) does not exist as 
a quantum theory. In the strict L = limit, and for massless fermions (or fermions with a finite mass), the theory 
does not have a ground state. In this paper, we study finite-L, N ^ oo limit of ([l]) which is a well-defined quantum 
theory. The Hermitian model p7l) will be useful as auxiliary to infer some general lessons about the finite-L case. 
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The /3 — )• oo limit of (j2ip is both intuitive and insightful. Since it is the limit of arbitrar- 
ily low temperatures, the fermionic boundary conditions should not matter, and indeed, this is 
transparent from (|2ip : the hyperbolic cosine increases unboundedly, and the trigonometric co- 
sine is bounded, hence the one-loop expression is dominated by the formerJl The ground-state 
energy of the system at one-loop order (or one-loop potential) can be deduced from the limit 
lim/3^oo 'S'ld.ioop//? = £^id,ioop[|^ab|] and the result is 

-Eld,l-loop[|^a6|] = 2^ - Xb\ - 4:Nf ^ \/|Xa -Xb|2 + m2, 
a<b a<b 

lim ^ld,l-loop[|Xab|] = (2 - mP) V \Xa - Xb\. (22) 
m— >0 ^ — ' 

a<b 

This is an intuitive and simple resul10 and various remarks are in order regarding the chiral 
(m = 0) limit on R. 

1. For A'^^ > 1/2, £'id,i-ioop[|2;a6|] is unbounded from below and eigenvalues mutually repel each 
other. The classical minima corresponding to the space of commuting triples = 
are unstable against perturbative quantum fluctuations. This means, at the L = limit, 
the Hermitian matrix quantum mechanics ()17p does not have a ground state. 

2. At finite L, since the target space of eigenvalues is compact three-torus with size 1/L, the 
repulsion implies that eigenvalues will uniformly distribute over T^. This implies unbroken 
center symmetry in the N = oo limit and the theory obeys volume independence. 

3. For < iV^ < 1/2, E^d, i-ioop[|2^afe|] is bounded from below. The minimum is at \xab\ = 0. 
However, in this domain, one-loop analysis is not reliable, and there are quantum fluctuations 
of order X^jj^. At A'^^ = 1/2, m = 0, the ground-state energy is zero to all orders in 
perturbation theory due to supersymmetry. 

In the following, we study in more detail the phase structure of the four-dimensional theory by 
using (|2ip and estimates of the nonperturbative quantum fluctuations. We first explain Nj^ = 
(bosonic) and A'^^ = 1/2 (or one Majorana fermion) cases. We then study the case with A'^^ > 1, 
our main interest in this paper. 

®Also note that at /3 — s- oo limit, At can be gauged away. Hence, it should not appear in one-loop potential. 
^'^ We could have guessed this result by physical reasoning as follows: In the weak-coupling domain where 
Aid ~ Aid/l^abP ^ 1 is small and one-loop analysis is reliable, the Lagrangian (|17|l reduces to a collection of 
bosonic and fermionic harmonic oscillators: 

L ~ i|9t^"''l^ + ^\^a — xtf\X^''\^ + fermionic oscillators . (23) 

The eigenvalue differences of background X matrices are identified with the frequencies of harmonic oscillators, 
= — Xb\ and oifj = \/\xa — Xb\'^ + TV?. In the chiral limit, cjf^ = uJab- In gauge quantum mechanics, there 
are 2 X]a<i) rnany massive bosonic fluctuations and 2 x iNf X]a<6 niany massive fermionic fluctuations. Hence, the 
background dependence of the ground-state energy of the system is 

E= k"^'- i'^'' = £^ld,l-loop[15a(,|] , (24) 

bosons fcrmions 

which is given in (|22|) . 
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3.2.1 Nj? = (bosonic) 

The one-loop action (j2ip generates an attraction between eigenvalues at long distance and 
the center symmetry is broken. More precisely, there are A^^ saddles related to each other by 
center conjugations. The tunneling between these saddles is suppressed in the large-iV limit 
and the theory is in a center-broken phase. However the eigenvalues do not collapse to a single 
point due to the nonperturbative quantum fluctuations; the clump has a finite size of order 
A}^^ = [A4rf(l/L)]V3/L which is 0(iV°) in the large-iV limit. 

3.2.2 Nj? = 1/2 (single Majorana) 

The super symmetric theory {Nj^ = 1/2 and m = 0) has been studied extensively, because its 
maximally supersymmetric cousin has a dual description as the DO-brane system. In this case the 
one-loop potential falls off exponentially, 

Sld,l-loop^^^^Pi-l^\^a- Xb\) , (25) 
a<b 

and the ground-state energy is -Eid.i-ioopll^abl] = li™/3-s>oo 'S'ld i_ioop//^ ~^ 0. This follows from 
the cancellation of the attractions and repulsions between eigenvalues due to supersymmetry. 
Consequently, there exists a space of flat directions. Therefore, once eigenvalues are well separated, 
they will propagate freely like the gas of DO-branes to all order in perturbation theory. 

In order to study nonperturbative aspects, it is more useful in this case to recall that the theory 
is formulated on X M, and to discuss the system in Hamiltonian formulation. As asserted above, 
to all orders in perturbation theory, the theory has a moduli space, and all possible realizations 
of center symmetry are possible. However, non-perturbatively, this is not the case. Quantization 
of the zero-momentum bosonic modes gives rise to a discrete spectrum and a gap. 

Witten studied the gauge quantum mechanics for finite- within Born-Oppenheimer approx- 
imation and showed that the bosonic ground-state wave function (ignoring fermionic zero modes 
which are not crucial in what follows) is constant 

^o{ei,...,eN-i) = i (26) 

and the excited states have an energy gap [3S]. The center symmetry and volume independence 
was not discussed in Ref.[38]; however, its results have natural implications which apply to our 
discussion. In particular, we can introduce an eigenvalue distribution function p{6) measuring the 
density of eigenvalues. The density is everywhere non-negative and obeys / d?9p{6) = 1. Since 
the ground-state wave function spreads uniformly over the perturbative flat directions, 

and the center is unbroken. Unlike the discussion in Sec l3.2.1"l it is also unbroken at large-A^. 
Nonperturbatively, we have a unique saddle singlet under center conjugations, and this is the 
main difference with respect to purely bosonic theory which has A^'^ saddles. 

When the fermion is massive, there is a subtlety due to order of limits of small mass vs. small-L 
analogous to the discussion in Sec. 13.1.21 with similar conclusions. 

Finally, at large- A^ , this theory has a metastable bound state of eigenvalues [l3] (with diverging 
lifetime as A^ — ?• oo) which is analogous to the one in maximally supersymmetric theory [46] 
corresponding to the black zero-brane in type IIA supergravity. 
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3.2.3 Nf > 1 

First consider the case with m = 0. From the one-loop effective action (j2ip . it is apparent 
that repulsive force coming from fermions dominates at long distance. In the Hmit where is 
shrunk to zero size, the resulting Hermitian quantum matrix model is not well defined; the one- 
loop potential is unbounded from below for large eigenvalue separations. For finite x M, the 
eigenvalues can no longer run off to infinity; instead they spread over the dual uniformly. We 
expect the center-symmetric phase to continue upon compactification, down to x Sj^ so long 

1/3 

as A^^ 13 ^ 1. In this domain, the center symmetry is intact, and in the large- limit, volume 
independence must hold. 

1/3 

In the high temperature limit, A^^^ /3 ^ 1, fermions decouple due to thermal mass, eigenvalues 
clump, and a metastable bound state of eigenvalues exists. Let us assume |/3Ax|, |/3Aa| <^ 1, 
where Ax = ma,Xa^b{\xa — Xb\} and Aa = niaxa^b{\(^a — Ofol}- Then the second term in the right 
hand side of (|2ip is negligible, due to the large thermal mass of fermions, and the effective action 
becomes that of the zero-dimensional bosonic matrix model. 



where we used the identification x'^ = Oa- This potential produces attraction between eigenvalues. 
This model is studied extensively, and taking into account nonperturbative effects, a bound state 
of eigenvalues exists and it satisfies conditions |/3Aa;|, |/3Aa| <^ 1. In this domain, (Zat)^ center 
symmetry is completely broken. 

Generalization to nonzero m is straightforward and follows from the discussion of Sec l3.1.3] 
on T^. At sufficiently small 0{N^) volume such that m <^ A|^^^, the effect of mass is small, and 
hence attraction in one-loop action is overwhelmed by nonperturbative quantum fluctuations. In 
this domain, the (^n)^ center symmetry is intact. 

4 Lattice model and Monte-Carlo simulation 

In Sec J3.H we explained on continuum and x why the Eguchi-Kawai reduction holds 
for the QCD(Adj) at zero and finite temperature, based on perturbative-loop analysis, supple- 
mented crucially with the estimates of nonperturbative quantum fluctuations. Below, we study 
the unitary matrix model and one-dimensional lattice model by using Monte Carlo simulation. 

It is hard to implement the thermal model on a computer since the temporal direction is not 
reduced. Moreover, in order to describe phenomena typical to the finite-temperature system, we 
need to take the effective spatial volum^^ sufficiently large compared to iVj. As for a lattice 
model with an isotropic lattice spacing, roughly speaking, we need to take the spatial lattice size 

The space size in its ordinary sense is fixed. By developing perturbation theory around a center symmetric 
background, say for simpUcity with only one-dimension compact, we observe that both the lowest states and the 
spacing between the states is suppressed by a factor of A'' and is given by 2n /{LN) as opposed to the usual (center- 
broken) KK-spectrum where the level spacing is 27r/L . In other words, to all orders in perturbation theory, the 
effective space size is enhanced into Leff = LN. (See the discussion in [^Hi for QCD(Adj) and related discussions in 
a review [36j for the quenched Eguchi-Kawai model (QEK) and the twisted Eguchi-Kawai model (TEK).) This is 
the perturbative essence of volume independence, and in the sense of neutral sector observables, decompactification 
limit can be reached at A*' — >■ cxd while keeping L fixed. 



S- 



Id, 1— loop 





a<b 
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Ns to be twice larger than the temporal one A'^^, Ns ^ 2Nt. Since the A'^^'^ is related to the 
matrix size N for the large- A'^ reduced models in the sense described in footnote [TTl we need to 
take large enough to satisfy this condition. When we use the original Eguchi-Kawai reduction 
for the spatial directions, finite- A^ correction behaves as 1/A^ rather than as we will show 

numerically. This fact indicates that the effective lattice size scales as A'^^'^ ~ N'^^'^, and it is not 
practical for numerical simulations. Therefore, to gain more effective spatial volume, we impose a 
twisted boundary condition on the spatial directions. For the latter, the effective volume is more 
enhanced, and volume independent domain can be reached more quickly. We call the QCD(Adj) 
with the twisted boundary condition as TQCD(Adj). The introduction of TQCD(Adj), which is 
algorithmically more convenient, is our main improvement over the one-site model of Ref . |32j . 

To explain the efficiency of the TQCD(Adj), we start with the zero temperature case, that is, 
the single-site model with periodic boundary condition along the temporal direction on fermions. 
Then we apply the twist to the finite-temperature case, that is, the one-dimensional lattice model 
with an antiperiodic boundary condition on fermions along the temporal direction, and then see 
the Eguchi-Kawai reduction holds in this case. Especially we show this reduced model at finite- 
temperature describes the confinement /deconfinement transition. 

4.1 Single-site theories 

4.1.1 QCD(Adj) on 1^ lattice 

The single-site matrix model can be obtained from a four-dimensional lattice gauge theory by 
reducing the number of lattice sites to one in all directions [32] . The action is 



Sod = -2bNRe Tr I V^V^^V^Vj J + Sp, 



(29) 



where (^u = 1,2,3,4) corresponds to the link variable in the four-dimensional theory. The 
inverse 't Hooft coupling constant h should be chosen appropriately depending on the lattice 
spacing a. The fermionic part Sp is obtained as dimensional reduction of the Wilson-Dirac 
fermion term 



= ^ vi/V'/ - K J] {^5/(1 - 7m)^m^/^/I + ^^/(i + i^)yl^fy^ ■ (30) 



f=i \ i=i 

The hopping parameter k can be expressed as 



K= , (31) 

8 + 2amo ' ^ ' 



where niQ is the bare mass. 

This action has a {"^n)^ center symmetry 

^ e^™''/^^^ (n^ = 0, 1, • • • , AT - 1). (32) 

If this symmetry is not broken, then the model is equivalent to the translationally invariant 
subsector of lattice theory with an arbitrary number of sites, including an infinite lattice limit. 
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Although detailed analytic evaluation of one-loop effective potential depends on the choice of 
lattice fermions0 intuitively, the absence of the center-symmetry breaking phase follows closely 
the discussion on continuum T^. The discussion on can easily be generalized to the 1^ lattice. 
The role of the compactification scale L is replaced by the lattice spacing a in the one-site model. 
The main lesson that we learn is that the center symmetry on the 1^ lattice model is in fact much 
more robust than the center on continuum T^. 

The one-loop action for the one-site theory in the classical background of commuting Wilson 



lines [V^, V^] = where V"^ = diag(e* 



, e* ) is given by 



Si- 



loop 



2 E t (f ) ] - ^^F E U t < + ( 

a<b /i=l a<b \ M=l 



oab , 



The first term is induced by gauge fluctuations and leads to eigenvalue attraction [13]. Geomet- 
rically, 



V 
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a 



je?, 



(34) 



is the separation between two eigenvalues of the Wilson line in the fj, direction. = X]^=i(^^'')^ 
is the spectrum of massive gauge fluctuations {W bosons), familiar from the usual D-brane pictures 
as the spectrum of open strings ending on branes, where eigenvalues of Wilson line are identified 



with branes. The second term, proportional to NP , is induced by fermionic fluctuations, and is 



equal to -4iVf Ea<6 log 



where M'j[9'^^] is the spectrum of fermions. 



Eq. (j33p may be rewritten in a form mimicking the continuum expression on given in ([T 



S 



\nab 
1-loop [C^ 



]=2E log 



a<b 



-mf log 

a<b 



4 2('pab\2cpab\2 

E«)' + E ' 



(1=1 



2(1 -|- moa) 



+ 



(1 -I- moa) 
(35) 



where we have subtracted a holonomy-independent constant term. This expression can be Poisson 
resummed and be written in terms of Wilson lines as 



Si-ioop[0l!'] = E E 

a<b neZ4\{0} 



''P^(moa) , 



(36) 



as discussed in the Appendix [XI The only difference with respect to the continuum expression on 
is that 

Pi^'*^(0) = e„^pf (0), (37) 

where is an enhancement factor of a one-site model over continuum T^, for both bosonic and 
fermionic contribution. In general, due to peculiarities of the dispersion relation of the one-site 
model, the center stability is further enhanced on the one-site model with respect to continuum 
T^. Asymptotically, for |n| » 1, Cft — > 1 as expected on physical grounds, by just inspecting the 
dispersion relations. 
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Here, we use Dirac- Wilson fermion with Wilson parameter r = 1, [33] uses overlap fermions, and also discusses 
naive fermions. 
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The fermionic contribution to is numerically sizeable and has interesting implications. For 
example, on T^, start with mL = oo. Following the discussion of the Appendix [XJ the singly- 
winding Wilson lines do get stabilized (in perturbation theory) at mL = 2.027, whereas for the 
one-site model, the same phenomena take place at moo = 9.3. This is due to the fact that in the 
domain of heavy fermion bare mass, the last term in ()35p . which may roughly be viewed as an 
"effective mass" m'^g, is suppressed with respect to the bare mass, meff ~ y/mo/a. For A^^ = 1 
theory, the fermions' contribution dominates, leading to eigenvalue repulsion, and unbroken center 
symmetry for the theory defined on a 1^ lattice, or EK reduction of QCD(Adj) to a single-site 
lattice. 

If one-loop perturbation theory was the whole description, this would be the transition to a 
(Z2)^ restored phase. However, as explained above, many phases can coexist and it is beyond 
perturbation to determine which phase is chosen in the end. Interestingly, the point where the 
center symmetry (partially) restores in simulation is very close to the prediction by perturbation 
theory k ~ 0.037. This is approximately the value of k where Ref.[32] observed the full center 
restoration, not just Z2. In Fig.fH we can also see the restoration around the same value. Although 
we could not see partial breaking clearly by measuring the observable employed in Ref.[32], the 
discrepancy between in large-K start and the one in small-K start may indicate the existence 

of a partial breaking phase. 

4.1.2 TQCD(Adj) at zero temperature 

In compactified QCD(Adj), as we will see explicitly from simulations, the finite- corrections 
turn out to be order 1/A^, as opposed to the perturbative expectation on R'* [3l]. As explained in 
the results section of Sec Jl.ll there are two plausibly related explanations for this behavior. One is 
related to the discussion of effective volume in the reduced model. In the reduced model, A^ serves 
the role of an emergent spacetime volume, at least in a perturbative description in finite volume 
around a center-symmetric configuration. Finite- A^ corrections should scale as finite-volume cor- 
rections. However, what is not always clear is the factor A^^ via which volume enhancement takes 
place Ves ~ NPV, El and p may in fact be determined nonperturbatively. 

In compact space, one cannot gauge away zero-momentum modes, and these modes are crucial 
in studying perturbation theory in finite volume. In perturbation theory, the spectrum of the 
theory relies on the background for the Wilson lines. If theory has massless adjoint fermions, 
there will also be fermionic zero modes in the spectrum. Typically, there are order A^ light or 
massless bosonic and order A^ fermionic zero modes, which may generate nonperturbative 1/A^ 
effects [35] . 

Both problems can simultaneously be solved and 1/A^ corrections can systematically be im- 
proved by using boundary conditions which cannot be obeyed by either bosonic or fermionic (if 
there are any) zero modes. This can be done by using the twisted boundary conditions of 't Hooft 
|47| . This idea is, of course, not new, and is used by Witten in Ref.[38] to lift the zero modes 
in M = 1 super Yang-Mills theory in the context of supersymmetric theories on X M, and by 
Gonzalez- Arroyo and Okawa jl5j in the context of large- A^ reduced models. 

^■^Of course, for QEK-like configuration, p = I and for TEK-like configurations, p = 2, see, for example, [36j . 
However, tliese deductions are in perturbation tlieory around particular backgrounds, and tlie determination of p 
is likely non-perturbative. 
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Our main observation can be summarized by using the fohowing pedagogical exercise. (The 
generalization to the theories that we use in simulations is straightforward. The prescription given 
below works equally well on lattice and continuum.) Let ^(xi, 2:2, X3, X4) denote either a unitary 
gauge field or an adjoint fermion field. We impose the following generalized boundary conditions 
on fields 



$(..., + L, ...) = B^^{. ..,x^,.. .)Bl . 



(38) 



For the particular case of one-site matrix models, we can set <I>(xi, X2, X3, X4) = ^ =constant. For 
the theory on the 1^ lattice, we consider two choices for 



pbc : 
Twist : 



Bi = C, 



N ' 



5 



Iv^^^v^ (39) 



where C /j^ and S are y N xy N (noncommuting) clock and shift matrices obeying C ^jjS, 



e ^ S /TfC /jf. A particular representation is 
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(40) 



where uj = exp(27ri/viV). 

The first case in ()39p is the original QCD(Adj) with periodic boundary conditions and does 
not lift any zero or light modes associated with holonomy, or fermions. In this case, finite-iV 
corrections turn out to be largest, of order 1/N. 

The second case in (j39p is a twist of QCD(Adj). The twist lifts all possible zero or light modes 
from the spectrum. In this case, finite- corrections turns out to be of order 1/N'^. 

The action of the theory with twisted boundary conditions can be turned into a theory with 
periodic boundary conditions and an action with an insertion of 't Hooft flux. This is our definition 
of "twisted" QCD(Adj) [or TQCD(Adj)]: 



Sod = -2bNRe Tr ( ^ VF^KF^yJ j + gp, 



(41) 



where Z^y is the twist factor. Geometrically, Z^^, is associated with the 't Hooft flux passing 
through the (/uz/) plaquette. Here we adopt "symmetric twist" 



Z, 



(42) 



As the area enclosed by fermionic "plaquette" terms is zero, the fiux passing through it is zero. 
Thus, fermionic action is unaltered. This procedure, apart from helping QCD(Adj) algorithmi- 
cally, also cures the global instability [171 \T9\ [T8] of the TEK model. We numerically compare the 
behaviors of finite- A^ corrections for QCD(Adj) and TQCD(Adj) below. 
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4.1.3 Numerical results for Nj' = 1 

We now discuss tlie Monte Carlo results for QCD(Adj) and TQCD(Adj) at zero temperature. 
We restrict our analysis to the case with a single Dirac fermion in adjoint representation!^ 

In Fig. [21 the expectation value for the absolute value of the Wilson loop (averaged over all 
directions) 

\W\^\Y.\^,1 (43) 

in the QCD(Adj) and the TQCD(Adj) at zero temperature is plottedQ For both the QCD(Adj) 
and TQCD(Adj), (|VK|) is of order 1/A^ and hence the (Z^v)^ symmetry is unbroken. (As already 
shown in [32], it is unbroken in a rather large parameter region.) The extent of the next-to- 
leading correction is not clear from this plot; we fit it by ~ c/N + d/N'^ for QCD(Adj) and 
{\W\) ~ c'/iV + d'/N^ for TQCD(Adj), where c, d, c' , d' are constants. 

In Fig. [31 expectation values of the plaquettes are plotted. From this plot, the finite- iV 
correction for the QCD(Adj) turns out to be of order 1/A^. On the other hand, the finite- 
correction for the TQCD(Adj) is of order as expected. 
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Figure 3: Expectation values of the plaque- 
tte in QCD(Adj) at 6 = 0.50, k = 0.09 and 
TQCD(Adj) at 6 = 0.50, k = 0.09, k = 0. 
The 1/A^ correction is of order 1/A^ for the 
QCD(Adj) and l/N^ for TQCD(Adj). 
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Figure 2: Expectation values of the Wilson loop 
in QCD(Adj) and TQCD(Adj) at 6 = 0.50, k = 
0.09 and k = (bosonic twisted Eguchi-Kawai 
model). Fitting curves are of the form c/N + 
d/N'^ for the former and c'/N + d' /N^ for the 
latter. TQCD(Adj) at 6 = 0.50, k = 0.09 and 
K = agree quite well, as expected because k = 
0.09, corresponds to a quite heavy fermion. 

^* We implemented the rational hybrid Monte Carlo (RHMC) algorithm [48] with the multi-mass conjugate 
gradient (CG) solver [49l. Numerical coefficients in the rational approximation necessary for the RHMC simulation 
was obtained by using the simulation code provided at [50| . 

^^We use absolute value of the Wilson line operator in small volume to distinguish a center-symmetric saddle 
point from a multi-saddle configurations for which (W) is non-vanishing at each saddle, but vanishes due to phase 
averaging over all saddles (which is permitted in quantum theory due to tunneling). Multi-saddle configurations, 
in the large-Af limit, lead to spontaneously breaking of the center-symmetry, whereas a center symmetric saddle 
continues to respect the center symmetry. 
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In Fig. m expectation values of Wilson lines near the phase transition are shown. The 

argument in Sec l4.l.l1 suggests the transition is of first order because several phases coexist. To 
confirm, we studied and observed hysteresis. We started simulation at k = 0.01 (small-K start) 
and n = 0.05 (large-K start), and gradually increased/decreased the value of k. At each point, we 
collected 500 - 2000 samples, which is enough to evaluate the expectation values. As can be seen 
from the plot, there is a clear hysteresis. Thus we conclude the transition is indeed of first order. 

We also studied the distribution of TriVi^yy) /N as in [32], and did not find partial breaking 
of the center symmetry. However, around k = 0.04 in Fig. 4, there is some discrepancy between 
in the large-K start and that in the small-K start. Hysteresis continues with a smaller gap 
within 0.0375 < n < 0.0475. This may be interpreted as a partial breakdown of center symmetry 
|51j , where the unbroken phase coexists with a partially broken phase. 
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Figure 4: Expectation values of the Wilson loop {\W\) in QCD(Adj) at 6 = 0.50, = 25. Clear 
hysteresis can be seen. 



4.2 One-dimensional lattice: Reduced model at finite temperature 

We first introduce a one-dimensional lattice formulation corresponding to the large- A^ reduced 
model at finite temperature and then apply two different types of twists to this model. An 
implication of volume independence is that one can study confinement/deconfinement transition 
on X on the equivalent unitary matrix model, corresponding to a 1^ x Nt lattice so long as 
{IinY center symmetry associated with the spatial cycles is not spontaneously broken. We indeed 
observe the confinement/deconfinement transition in the reduced model. For the comparision of 
1/A corrections, we also plot some numerical results for the QCD(Adj) (without twist) at finite 
temperature. 

4.2.1 (T)QCD(Adj) at finite temperature 

On the 1^ X Nt lattice or x continuum formulations, the relation between twisted boundary 
conditions, zero (or light) modes, and suppression of finite-A effects can systematically be studied. 

Let $(t, X, y, z) denote either a unitary gauge field or an adjoint fermion field. We impose, the 
following generalized boundary conditions on fields 

$(t, x + L,y,z) = A<^{t, X, y, z)A^ , 
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x,y + L,z) = B^{t, X, y, z)B^ , 
^t,x,y,z + L) = ^{t,x,y,z). (44) 

For the particular case of matrix models, we can set <^>(t, x, y, z) = ^{t). We consider three choices 
for A and B. 

pbc : A = 1n, B = lj\f 

Twist 1: A = C^<g)l^, B = S^'^1^ 

Twist 2 : A = Cn , B = Sn (45) 

where Cat and Sn are N x N clock and shift matrices defined earlier. 

The first case in (j45p is the original QCD(Adj) with periodic boundary conditions and does 
not lift any zero or light modes associated with holonomy, or fermions. In this case, finite-A^ 
corrections turn out to be of order 1 /N. 

The second case in (j45p is a partial twist of QCD(Adj). The twist lifts a fraction of 

bosonic light modes and fermionic zero modes (in cases where fermions are light). This can be 
seen by explicitly solving the boundary conditions. In this case, finite-A^ corrections turn out to 
be of order l/N^/"^ . This is explained below after introducing the twisted model. 

The third case in ()45p is a twist of QCD(Adj), which lifts all bosonic light modes and fermionic 
zero modes. This can be seen by explicitly solving the boundary conditions. In this case, finite- A^ 
corrections seem to be rather tame. 

We study the theory reduced to a one-dimensional SU(N) unitary matrix model 

Slat = -2bNY,Re Tr ^ Ut{t)Vi{t + l)u}{t)V^{t) + J] ZijV,{t)Vj{t)V^ {t)vj (t) j + Sp, (46) 
t \ i i<j J 

where Vi{t) {i = 1,2,3) corresponds to the link variable along the spatial direction in the four- 
dimensional theory. For the usual nontwisted model Zij = 1 and Za = 0, while, for twisted 
models, Zij is given bjj^fl 

Twist 1 : = = e^^*/^ {i < j), 

Twist 2 : = Z*i = e^^*/^ (^ < j)- (47) 

The number of sites At is related to the temperature Tby/3 = 1/T = aNt- The fermionic part 
Sp '^s given by 



Nf 



Sp = J2^fD\Pi;f, (48) 
/=i 

where Dw is the usual Wilson-Dirac operator. The explicit form of Sp after reducing the spatial 
directions is 

t ^ i=i 



^®In simulations, we use symmetric twist which is more efficient. The above boundary conditions can be modified 
to incorporate the symmetric twists. 
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-K - itwmit + i)ui{t) + + it)u}xt - i)i^{t - i)Ut{t - 1)} j • 

(49) 

We impose antiperiodic boundary condition for the fermions on the and generahzed bound- 
ary conditions given in (j45p on reduced directions. The action has a global center symmetry, which 
we split for convenience as {^nY ^ '^N-, 

Vi(t) ^ e2-/^y.(t), {Ut=i...Ut=N,)^e^^'"'{Ut=i...Ut=N,). (50) 

If the (Ztv)'^ symmetry is not spontaneously broken, then this model is equivalent to the one in 
the infinite spatial volume lattice. This implies the phase transition of infinite volume theory, 
associated with the realization of the temporal Z^v factor, can be studied by using the unitary 
matrix model. 

The relation between the twists, number of light modes, and the observed form of the finite- A'" 
corrections is 



Twist 


Zero modes 


finite— corr. 


none 


N 


1/iV 


e viv 


Vn 


l/Ar3/2 


e N 


none 





(51) 



There is a nice geometric interpretation for twist 1 in terms of a classical background, and 
the foliation of noncommutative plane. This gives, in perturbation theory, that effective volume 
should scale as y ~ N^/"^. However, for twist 2, there is no classical background solution. Because 
of strong quantum fluctuations, a classical background cannot be written. Of course, this is not 
a concern. 

The emergence case of twist 1 can be explained in perturbation theory as 

follows: For simplicity, let us again consider the twist Z12 = e^'^*/^^, Z13 = Z23 = 1. A natural 
candidate of the ground state is 

Vi = C^®l^, V2 = S^®1^, ^3 = 1^®^^. (52) 

This is because this configuration satisfies ViVj = Z*-VjVi, and at the same time eigenvalues 
spread as uniformly as possible. This configuration keeps the (Z subgroup of the center 
symmetry, which is enough for the Eguchi-Kawai reduction to work. Then, along Vi and V2 
directions ^/N sites arise as a fuzzy torus, similar to TEK, and \/iV sites emerges along the V3 
direction similar to QEK. So this configuration corresponds to V^^ ~ JSf^/'^ lattice sites. If one 
views finite- A^ corrections around N = 00 analogous of finite-volume corrections in compactified 
theories, then, it is expected that finite- A^ corrections in the one-dimensional reduced model should 
be a power series expansion m 

4.2.2 Monte Carlo simulation for A'^^ = 1 

Here we show the Monte Carlo results for the one-dimensional lattice model for A^^ = 1. The 
numerical algorithm adopted is the same as the one in the previous subsection. 
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First let us see that the (Ztv)^ center symmetry is not broken. In Fig.O we plot the expectation 
value of an averaged Wilson line given by 

1 3 Nt 

at b = 0.5 and various k, and A'^t = 1 for QCD(Adj). As expected, the is not broken when 

quarks are sufficiently light for both models. We can see similar behavior for the TQCD(Adj) 
except for the finite-A^ corrections. From Figs. [6] and [71 it is clearly seen that goes to zero 

as for the QCD(Adj) while l/iV3/4 f^^. ^.j^g TQCD(Adj). That the Wilson hue behaves as 
{\W\) ~ iV^3/4 ^ ^j^g TQCD(Adj) is the same as what happened in the model in the 

previous section. The same pattern can also be observed with A^j > 1 . 
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Figure 5: Expectation value Figure 6: at 6 = 0.50, k = Figure 7: at 6 = 0.50, k 

of the spatial Wilson line in 0.10, iV* = 1 for (nontwisted) 0.10, A'^t = 1 for TQCD(Adj). 

QCD(Adj) (without twist) at QCD(Adj). 

b = 0.50. When the quark is 

lighter, the Wilson line becomes 

zero at large- A?"; that is, (Zat)^ 

center symmetry is not broken. 

In Fig. [HI the plaquette is plotted for the QCD(Adj) and TQCD(Adj) with two types of twist. 
The Ansatz const. + const. /N^^"^ for the twist 1 is consistent with the data, and by assuming it, 
the extrapolated value at A^ = oo agrees with the one obtained from the nontwisted model. For 
twist 2, const. + const. /N"^ is consistent with the data. The similar behavior can be seen for the 
expectation value of the Polyakov loop defined by 



\P\ 



1 



Nt 

TrX{Ut{t) 

t=i 



(54) 



We observe that by using twist 2, we can suppress finite- A^ corrections more. We notice that, up to 
A^ = 25, we do not observe the jump in the expectation value of the plaquette, which corresponds 
to the bulk transition. 

Now, let us consider the confinement /deconfinement phase transition in the large- A^ reduced 
model at finite temperature. We take k = 0.10. From the experience in pure bosonic Yang-Mills 
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theory, which is studied extensively in [391 140j . it is known that one has to take the number of 
spatial links Ns to be large (roughly Ng > 2Nt) in order to see the deconfinement transition clearly. 
For QCD(Adj) without twist this is rather severe constraint, because A'^^''^ is related to the number 
of colors A^ by A^^'^''^ ~ A^^/^. As a resuh, at Nt = 2, we could not observe the transition below 
A^ = 30, although we could observe a clear transition for Nt = 1. In contrast, in TQCD(Adj), the 
transition can be seen at small values of A^. In Fig. [10] we plot the expectation value of the Wilson 
line and the Polyakov loop for At = 2. We can clearly see a jump of the expectation values of the 
Polyakov loop around b = 0.33 for Nt = 2, which corresponds to the confinement/deconfinement 
transition^ This result is consistent with the one given by large volume simulations for pure 
Yang-Mills theory (for example, SU{8) YM theory simulated on the 10^ x 5 lattice gives, in the 
extrapolated infinite volume limit, be ~ 0.34 [39[ I40|). confirming the finite-temperature version 
of equivalence nonperturbatively. 

Before closing this section, let us give an estimate for the the physical temperature of the 
transition. By using the 2-loop beta function of the bosonic Yang-Mills theory (neglecting fermions 
because they are heavy), b is related to be the temperature T and lattice lambda parameter A^at 
as 

_ J_ f ^ Y^/^^^ /127r2 
Alat~ Nt\247r^V '""H HA 

where A = 1/(26). If we take Alat to be of order 0{lMeV) as in the QCD with SU{3), by sub- 
stituting it to ([55]) . the transition temperature turns out to be of order O(lOOMeV), as expected. 

^^In order to take the continuum limit, we need to analyze large Nt. However, to do this, we also need to take A'" 
large such that Ns^^ > 2Nt. It is beyond our current numerical resources and we leave detailed analysis for larger 
Nt for future work. 
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5 Stabilizing noncommutative Yang-Mills theory 



There is a well-know perturbative equivalence between the TEK model and noncommutative 
YM theory. Nonperturbatively, this relation is problematic due to a global instability |17[ [TE[ [1^ 
[20] , In the perturbative description, the twist-eater background is used to generate a fuzzy torus, 
the noncommutative base space of target theory. The spontaneous center-symmetry breaking in 
the TEK model is associated with the spontaneous collapse of the noncommutative torus. This 
instability, in effect, is related to the center-symmetry breaking instability of the TEK model. For 
a nice discussion of the relation between noncommutative theories and matrix models, see, e.g., 

m- 

In this work, we suggest that TQCD(Adj) can be used to define Yang-Mills theory on non- 
commutative space (for recent developments, see e.g. [53]). 

As asserted above, TEK construction is problematic due to instability of the twist-eater back- 
ground \19\ [20} [T71 [T8] . Similar construction with the fuzzy sphere also fails |20j . In noncommu- 
tative field theory, this instability corresponds to tachyonic modes in gluon propagator \5A\ I55j. 

When adjoint fermions are introduced, the situation is different. The center symmetry is sta- 
bilized even at the one-site lattice and the twist-eater configuration is stable. Hence TQCD(Adj) 
with appropriate parameter scalings can serve as a nonperturbative formulation of noncommu- 
tative Yang-Mills theory with adjoint fermions. The interesting point is that the background is 
stable even if fermions are very heavy. In such a situation, low energy physics is almost a bosonic 
one, but instability in the one-site model is removed thanks to fermion induced stabilizing effects. 
We expect that the same effect can also be achieved by introducing a twist to the double-trace 
deformation of the Eguchi-Kawai model [31], which is purely bosonic. 

6 Conclusions and Discussions 

In this paper, we investigated the volume reduction for large-A^ gauge theory with adjoint 
fermions [7]- We used perturbative one- loop analysis crucially supplemented with the estimates 
of nonperturbative quantum fiuctuations to analytically explain the zero temperature result of 
Ref. [32]. We have shown that for heavy fermions there exists a large-small volume equivalence, 
with an intermediate volume dependent phase. In the sense of dynamics, the theory exhibits a 
4d-0d-4d cascade of dimensions. We have used the small volume phase to extract results for large 
volume theory. 

Next, we generalized the volume independence to finite temperature, and confirmed it numer- 
ically. In effect, we have constructed a one-dimensional lattice model and analyzed it by using 
Monte Carlo simulation to directly show that the center symmetry is preserved for a wide param- 
eter region including heavy fermions. In particular, we introduced a twisted version of the large- 
reduced theory with adjoint fermions for numerical efficiency and then succeeded in observing the 
confinement /deconfinement transition. The temperature agrees with large volume simulations of 
pure Yang-Mills theory. We also argued that TQCD(Adj) with appropriate parameter scaling can 
serve as a nonperturbative formulation of noncommutative Yang-Mills theory. 

There are several directions for future research. Our discussion of QCD(Adj) with the A^^^ = 1 
fermion can be generalized to address mass spectrum and other nonperturbative aspects. By 
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introducing more fermion flavors, our reduced models may be used to address the conformality 
or confinement problem, to address the determination of the lower boundary of the conformal 
window and perhaps to study models relevant to the technicolor scenario. 

In this work, we have seen that due to nonperturbative quantum fluctuations, the restoration 
of full center symmetry occurs not at niLN ~ few, but rriL ~ few. This suggests that the 
number of double-trace operators suggested in [31] for theories on T"* and x M is a conservative 
overestimation, and might be reduced considerably without spoiling unbroken center symmetry. 
If so, since deformed reduced theories are purely bosonic, the simulation cost becomes cheaper 
and it can be a practical tool. If it works, it can be used as a template to study QCD with 
fundamental matter. 
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A QCD(Adj) on continuum at finite-L and one- loop analysis 

In this appendix, we will rewrite the one-loop potential on small-T^ (|14p in terms of Wilson lines 
by using Poisson resummation. Poisson resummation is a duality which maps a sum over KK- 
momenta given in (|14p to a sum over a winding number. The two are equivalent expressions, and 
some aspect of physics is more transparent in one of the two. 

Since (|14|) is periodic under 9'^^ — ?• 6'^ + 27r, it can be Fourier expanded: 

Si^iooAOf] = Yl E ^''"■"PrtimL). (56) 

a<bneZ-i\{0} 

For a given winding number n = (ni, . . . , n^), 

a<b a,b o, 

It is also useful to express this sum in terms of trace over adjoint representation matrices, 
given by 

f^adj(n) = (Fi"^ . . . V2') ® (^r • • • K"')^ , (58) 



28 



or equivalently, 







'i-NxN 

























Clearly, trr2adj(^) = |tr (Vj"^ • • • V^^**)!^. Equation (jl4p can dually be rewritten as a sum over 
winding modes 

= ^ E <^ (^"' • • • ^4™^)!' - ^) ' (60) 

n6Z4\{0} 

where is interpreted as the mass square of the Wilson line with winding number n. Unlike 
the similar sums appearing in one-loop effective potential with 1 < d < 3 which are 

absolutely convergent [7], the series ()60p is conditionally converaen^^ This is, of course, physical 
and related to nontrivial infrared (IR) aspects of the theory which we discuss below. 



A.l IR singularities, conditional convergence, and noncommutative saddles 

The series in (jl4p and (j60p are equivalent expressions, related to each other via Poisson resum- 
mation. Consider the zero KK-momenta subsector of ()14p . Since eigenvalue separation has an 
interpretation as momentum, (jl4p exhibits an IR singularity whenever = for some 

a, 6. This IR problem is also manifest in lattice one-site one- loop action (j33p . The series (j60p is 
conditionally convergent, and whenever two eigenvalues are coincident, '}2n&z-^\{o} exhibits log- 
arithmic IR divergences due to modes with large-winding number |n| — >• oo. Winding modes 
have an interpretation in terms of spacetime distance and this is the same IR problem as in 
(jl4p . The physical interpretation of divergence is as follows: Whenever two (or more) eigenvalues 
are coincident, there are (in perturbation theory) massless modes (analog of W bosons or open 
strings). IR divergence comes about because we have integrated out these massless modes that 
we should have kept in the correct description of long-distance dynamics of the effective theory. 

Whenever two (or more) X^^=i(^)5l'')^ = 0, the zeroth order assumption that one can expand 
the fluctuations around commuting saddles (HD is incorrect . There are circumstances where 
commutative saddle points of the classical theory may be a good description at one-loop order 
in perturbation theory for some range of m. For example, if m = 0, then, the one-loop effective 
action (fH|) and ([60]) reduces to 

. ^ ^ , ^ (27rA;. + Off 



5i-ioop[<] ={2-ANf)Y^ E 



a<b fci ,...,fe4 \ At=l 



^*In (|60p . the subtraction of the constant divergent term is related to the absence of log divergent and holonomy 
independent X]a=6 terms in the first line, as well as (|14|) . With this subtraction, any singularity that may appear 
in the sum (|60p or its dual (|14p is physical IR-singularity. IR-aspects of one-loop effective action can be examined 
in either picture. 
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= {-l + 2Nf)^ J2 ^(trO,dj(n)-iV). (61) 

nGZ4\{0} ' ' 

This one-loop action is bounded from below and generates a repulsion between eigenvalues. 

If m = oo is sufficiently large, then the one-loop action ()14p and (j60p is unbounded from below 
and leads to an attraction between eigenvalues. A configuration where all eigenvalues clump and 
center is broken is the minimum. However, the action is IR divergent there. This is an artifact of 
perturbation theory. Whenever eigenvalues are close to each other, then the relevant scale in the 
theory is Aq^^, and the dynamics is strongly coupled for eigenvalues in the \xah\ ^ Aq^^ domain 
where one should not use perturbation theory. 

At one-loop order in perturbation theory, the actions (fH|) and (|60|) as well as the one-site 
version (j33|) realize all the saddles conjectured to exist in Ref . [56] . At m = 0, the leading fluc- 
tuations are quadratic or Gaussian; at m = oo, the leading fluctuations are quartic according 
to the classiflcation of Ref. [56]. As mL is dialed, all interesting saddles with varying number of 
quadratic and quartic fluctuations appear in massive QCD(Adj). In Refs.[5B] [T3]. only the two 
extreme cases were shown to exist in pure Yang-Mills theory in d dimensions, as d is varied. 

Since the one-loop action has IR singularities, the way to obtain the set of saddle points 
requires some care. 

1. Introduce an auxiliary IR cutoff ^ m, modifying gauge contribution in ([T3D to X]^=i(2'''"^a'+ 
6'^*)^L~^ -|- [35] • At any flnite but inflnitesimal value of /i/R, we can sensibly compare 
the one- loop effective action of different saddles, and flnd the global minima at a given value 

of mL. At m = oo, the global minima are at = 1 (and its center conjugates). This can 
be studied by using the Hermitian matrix model (|15p . 

2. Assume that the global minimum is a /c-bunch configuration of Wilson line phases. The 
nonperturbative width of each clump, due to quantum fiuctuations is determined by zero- 
dimensional matrix model. Interactions between different clumps are well approximated by 
one-loop effective action, and are repulsive. Dynamics inside each clump is approximated 
by the SU{N/k) matrix model, 

/ \ 

where we put primes in order to emphasize that matrices are {N/k) x (N/k). Then, 't Hooft 
coupling effectively becomes Xod/k, and hence the nonperturbative width of each clump is 
~ (Aod/fc)V4 . 

3. At relatively large values of fermion mass, mL ~ 1, one may expect that only a small sub- 
group of (Ztv)^ symmetry persists. However, if k is sufficiently large that distance between 
bunches becomes smaller than the fluctuation scale, 27r/(Lfc^/^) < Aq^^ = X\^/L, interaction 
between nearby bunches cannot be evaluated by one-loop approximation, and distinction of 
bunches becomes obscure. This implies that the /c-bunch phase is indistinguishable from the 
uniform phase. (See Fig JTT]) . 
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bunched phase (k=6^) uniform phase 




Figure 11: If k is large enough so that the distance between two nearest bunches 2Tr/{Lk^^^) 
becomes smaller than the quantum fluctuation scale Aq^^, interaction between these two bunches 
cannot be evaluated by using one-loop effective action. Quantum fluctuation turns fc-bunch phase 
into the uniform phase. 
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